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Preface 


In the late 1990s research was being performed at the Kennedy Space Center to develop in situ 
resource utilization technology for Mars. One study topic was the generation and storage of 
liquid oxygen (LOX) obtained from the atmosphere or regolith, but the transfer of this 
commodity was of concern. Mechanical LOX pumps were deemed potentially too heavy and 
unreliable for an autonomous mission to Mars, and alternatives were sought. One option was to 
use the paramagnetic property of LOX, which is significant enough that small quantities of LOX 
can be suspended against earth gravity with a rare earth magnet. With this application in mind, a 
small, internally funded project was initiated at the Kennedy Space Center in 2000 to study the 
use of pulsed magnetic fields to pump LOX. 

Proof-of-concept testing verified the LOX pumping predictions and resulted in a journal 
publication [1]. Numerous small coils were fabricated on cryogenic flow lines and used to 
produce intense, short-duration magnetic fields resulting in dramatic motion of the LOX. In 
addition, effort was expended on modeling the paramagnetic forces in the LOX, which required 
modeling the magnetic field generated by the coils and the coil inductance, allowing the current 
versus time to be predicted and compared to experiment. 

During the coil inductance modeling, we generated a method by which the inductance of an air- 
core solenoid could be computed efficiently. We wrote a short manuscript documenting this 
work, but did not proceed with publication. Instead, we chose to release our software (primarily 
generated by one of the authors, John Lane) with the manuscript attached as a reference 
document. The software and manuscript were made public through the Open Channel 
Foundation, an organization that publishes software from academic and research institutions. The 
link for this site is given below. 

http://www.openchannelfoundation.org/proiect/view docs.nhn? group id=288 

Through this link, our manuscript has been downloaded multiple times, has been used by several 
researchers, but to our knowledge, has not been cited. We suspect that this may be due to the 
lack of a clear citation and have decided that one should be provided. We consider it too late to 
seek peer-reviewed publication of this work and instead have decided to issue a NASA Technical 
Memorandum, with the same title and authors as the original manuscript. This will allow our 
work to be located by the search engines and will provide a reliable, long-term source and 
citation. The original manuscript, reformatted, but not in any other way altered, appears on the 
subsequent pages of this memorandum. 


[ 1 ] Youngquist, Robert C., Immer, Christopher D., Lane, John E., and Simpson, James C., “Dynamics of a Finite 
Liquid Oxygen Column in a Pulsed Magnetic Field,” IEEE Transactions on Magnetics, 39(4), July 4, 2003, 
pp. 2068-2073. 
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1 INTRODUCTION 

The pumping of liquid oxygen (LOX) by magnetic fields (B field), using an array of 
electromagnets, is a current topic of research and development at Kennedy Space Center, FL. 
Oxygen is paramagnetic so that LOX, like a ferrofluid, can be forced in the direction of a B field 
gradient. It is well known that liquid oxygen has a sufficient magnetic susceptibility that a strong 
magnetic gradient can lift it in the earth’s gravitational field. It has been proposed that this 
phenomenon can be utilized in transporting (i.e., pumping) LOX not only on earth, but on Mars 
and in the weightlessness of space. 

In order to design and evaluate such a magnetic pumping system, it is essential to compute the 
magnetic and force fields, as well as inductance, of various types of electromagnets (solenoids). 
In this application, it is assumed that the solenoids are air wrapped, and that the current is 
essentially time independent, 

2 B FIELD EVALUATION 

The magnetic field B due to a single current loop of radius a, as shown in Figure 1, can be 
computed by evaluating the curl of the magnetic vector potential A. With the current loop 
located in and symmetric with respect to the x-y plane, A is given by [1]: 

A (r 0 \ = Mo?L f 2g cos '<p' dqf (1) 

9 4/r * 0 +r 2 -2arsin<?cosp' 

where r, 6, and <j>' are spherical coordinates, i is the loop current, and po is the permeability of 
free space. The curl of Equation (1) can be simplified in cylindrical coordinates [2], or in 
Cartesian coordinates as follows [3]: 

^(w)=^ 2 ^( ( a 2 +R 2 )E(k 2 )-a 2 K(k 2 )) (2a) 

By(x,y,z)=^^ T {(a 2 +R 2 ) E (k 2 )-a 1 K(k 2 )) (2b) 

B A x ,y,z)=^{(a 2 -R 2 )E(k 2 )+a 1 K(k 2 )) (2c) 

where p 2 = X* +y 2 , R 2 =p 2 + z 2 , a 2 = a 2 + R 2 - lap, fi 2 = a 2 + R 2 + lap , A 2 = 1 - a 2 //? 2 , and 
C = p o i/tt. E(l?) and Kiji 2 ) are the complete elliptic integrals of the first and second kind, 
respectively. 
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z 



3 FORCE FIELD EVALUATION 


In order to estimate the force on a paramagnetic volume of LOX, it is necessary to evaluate the 
force equation [4]: 


F 



Mo 


K l+ Z. 


(B V) B 


(3) 


where again /Jo is the permeability of free space and X is the magnetic susceptibility of the 
medium. Equation (3) can be rewritten in matrix form as: 


where 
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It can be shown that D can be written as: 


D = PE(fc 2 )+Qar 2 K(& 2 ) 
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where P and Q are symmetric matrices. The five independent components of P and Q are given 
in Appendix A and Appendix B, in terms of thex,y, z field coordinates and loop radius a , where 

■?=*-?[ 3]- 

4 INDUCTANCE FORMULAS 

A common method of computing inductance of a solenoid in many engineering applications is to 
identify a handbook formula that most closely resembles the particular solenoid geometry at 
hand. Afterwards, a fine-tuning can be performed when building the coil by adding or 
subtracting loops of wire. 

For example, an estimate of the total inductance L of a single circular current loop is given by: 

L = o// 0 [ln(8fl/6)-|] (8) 

where a is the loop radius and b the wire cross-sectional radius, as shown in Figure 2. This 
estimate of self-inductance is good as long as bla « 1 . 



4.1 Stacked Single Layer 

Figure 3 diagrams a stack of M current loops, as would be found in a single layered solenoid. 

The handbook formula [5] for approximating inductance of this configuration is: 

(9) 


^/ 0 u 2 M 2 

ka 2 +(M -\) 2 d 2 


3 
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z 



Figure 3. M stacked current loops. 

4.2 Stacked Multiple Layers 

Figure 4 depicts a stack of M z current loops, configured as M x layers, as would be found in a 
multi-layered solenoid. The handbook formula [6] for approximating inductance of this 
configuration is: 


r 0 /+0.9 tg +0.32 / t+0.84/j)t 


(10a) 


l = {M z -\)d z t = 2M x b 
r o =o+(M x —\) b 


(10b) 
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4.3 Inductance From the Stored Energy 

Inductance L can be derived from the energy stored in a magnetic field due to a current carrying 
wire. In this case, L is equal to the integral of the product of current density J and magnetic 
vector potential A, divided by current squared [1 ]: 


L = 



01 ) 


For a single wire loop of wire diameter b and loop diameter a, the magnetic vector potential 
inside the wire can be approximated by [1]: 


A 


/V 

4n 


[(.- 4 1 

+ 2 

f'"- - 2)1 



b )\ 


( 12 ) 


with the assumption that bla « 1 . Assuming also that the current density is uniform (no high 
frequency skin effects, etc.), then J -i/Tcb 2 e . Integrating JA over the wire cross-section and 
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total length, as described by Equation (11), and dividing by i 2 leads to an approximation of self- 
inductance given by Equation (8). Note that due to symmetry, the volume integral becomes a line 

integral, ->4;r 2 aJ*p dp. 

Equation (8) describes an approximation to the total inductance of a single circular wire loop, 
which becomes exact as bla approaches zero. When multiple loops are considered, as would be 
found in a solenoid, the total inductance is the sum of two parts, self-inductance and mutual- 
inductance. Using Figure 4 as a general solenoid model, composed of M z stacks and M x layers, 
such that the total number of loops is equal to M x M z . 

The total solenoid inductance can be computed by summing the self-inductance and mutual 
inductance contributions of each m th , n [b loop (refer to Figure 5): 


( A /_ — 1 )/2 My - 1 

L = £ £ {L m „+M mn ) 

m=-(M : -l)/2 n = 0 

th th 

where the self-inductance of the m ,n loop is approximated by Equation (1 ): 

Kn=M o a „[H Sa n /b )-i] 

where a is replaced with a n = a + nd x . 


(13) 


(14) 


Mx 

A 



w=0 n=\ n=2 n= 2 
k = 0 k= 1 k=2 k= 2 


Figure 5. Solenoid model for general inductance calculation. 

The mutual inductance component M mn of each wire loop is determined by considering the 
magnetic vector potential of a single loop centered at the origin, in spherical coordinates: 

a _ AA' 

Wfl 2 + r 2 + 2<7rsin# 


(2-<*>W)-2£(l) 

<t> 


(15a) 


0 = 


Aar sin 6 


a 2 +r 2 + 2arsin0 


(15b) 
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where K(<p) and E(<p) are again the complete elliptic integrals of the first and second kind, 
respectively. 

Equations ( 1 5) can be used to describe the magnetic vector potential from each m , n wire loop, 
simply by translating the coordinate system. The r coordinate in Equations (15) describes the 
field at the m th , n ih wire loop due to the / h , k th wire loop, excluding the self-potential where j = m 
and k = n. To quantify this relationship, it can be shown that by following definitions given 
Figure 5, r in Equations (15) becomes: 


r jhm = 4 a l + ( m ~ j) 2 d; 

leading to the associated modification of Equation (15b): 


(P jkntn 


A 


K +a„Y +(m-j)-d: 


(16a) 


(16b) 


where a k = a + kd x . The magnetic vector potential from Equation (8a) for each m th , n lh wire loop 
now becomes: 


(^f>) jkmn 


Ao'V K + q n ) 2 + ( w - j ) 2 

K a„ 


((2 _ <Pjkmn )K( < Pjkmn ) “ ^^jkmrS) 


07) 


Finally, to express the mutual inductance of the m lh , n h wire loop, the integral of Equation (11) 
can be approximated by a simple product of the wire cross-section with a summation of 
Equation (17) over j and k, excluding j = m and k = n: 

.. (A/.-D/2 M t -I 

= = X I M +«.) 2 +(».-jf d- 08 ) 

^ j=-{M.-\)/2 A=0 
(excludingy=m, k-n) 


5 EXPERIMENTAL RESULTS AND DISCUSSION 

Magnetic field and force plots using two example coils are generated by employing 
Equations (2) and (4). The first coil has dimensions of a = 20 mm, b = 0.5 mm, d x = d z = 1 mm, 
M x = 20, and M z = 31. The B and force field plots are shown in Figures 6a and 6b. Figures 6c and 
6d are the B and force field plots of coil #7 in Table 1 . 
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F : (x, 0,z) 



F x (x, 0^) 



|F(*,0,z)| 

Figure 6b. Solenoid F-field plots from Equation (4) with dimensions a = 
20 mm, b = 0.5 mm, d x = d z = 1 mm, A7 V = 20, and M : = 31. 
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Table 1 . Dimensions of test coils [mm]. 


Coil# 

a 

b 

d x 

d : 

K 

M-_ 

1 

6.41 

0.051 

0.011 

0.013 

4 

150 

2 

4.57 

0.127 

0.240 

0.277 

27 

11 

3 

4.70 

0.255 

0.470 

0.540 

21 

33 

4 

4.70 

0.255 

0.540 

0.540 

27 

15 

5 

4.60 

0.140 

0.267 

0.280 

31 

19 

6 

1.67 

0.080 

0.170 

0.178 

17 

19 

7 

32.1 

0.362 

0.683 

0.757 

8 

83/84 


Equations (13), (14), and (1 8) make up a complete description of the total inductance of a 
solenoid described by Figure 5. In order to get some qualitative insight into the relative accuracy 
of this description of inductance, it is helpful to plot the self-inductance and mutual-inductance 
contributions, as shown in Figures 7 and 8. For these particular examples, the mutual inductance 
contribution quickly dominates as the number of turns increases. This implies that error resulting 
from the simplifying assumption of bla « 1 , used to compute the self-inductance term in 
Equation (14), will in most cases not be the primary source of error. Therefore, it can be 
concluded that accuracy of this method for computing solenoid inductance is dominated by the 
error introduced by allowing the integral in Equation (1 1) to be calculated as a product of the 
magnetic potential at the center of the wire cross-section times the cross-sectional area. 



Figure 7. Self- and mutual-inductance vs. M z for M x = 1 . 
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Figure 8. Self- and mutual-inductance vs. M x for M z = 1 . 

Figure 9 compares the calculated to the measured inductance of seven coils with dimensions 
given in Table 1 . “L Model #2” is equivalent to Equation (9); “L Model #3” is described by 
Equations (10); and “Z, Model #4” is represented by Equations (13), (14), and (18). Note that the 
simple approximation from L Model #3 seems to consistently track the more accurate calculation 
of L Model #4. Figure 1 0 is similar to Figure 9, but was generated using a Monte Carlo 
approach, comparing all inductance formula calculations to Model #4 results. 
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Figure 9. Measured vs. calculated inductance. 
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Figure 10. Monte Carlos test of inductance formulas. 


Even though these computational methods may seem “brute force” by traditional standards, the 
computational power of today’s desktop computers makes this approach very practical. 
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APPENDIX A. COMPONENTS OF MATRIX P 


[aV(3 Z W) + / >W-/)) 

-a 2 [p 4 {5x 2 +_y 2 )-2/? 2 z 2 (2.x 2 +y 2 )+3z 4 / 2 ') 
-R 4 ( 2x‘ 4 + 'fiy 1 +z 2 ))J 

|> 4 (3,> 2 -2z 2 )-S 4 (2S 2 + p 2 ) 
-2a 6 -2a 2 (2//-pV + 3z 4 )] 


W7 

+2z 2 (a 4 - 6a 2 /? 2 + /? 4 ) + z 4 (a 2 + /? 2 )] 


2^7 [aY+^jV-pV-*/)) 
+a 2 ( 2p 2 z 2 (x 2 + 2 y 2 ) - p 4 (x 2 + 5y 2 ) + 3 'f'z 4 ) 
+R 4 (f(x 2 +z 2 )- 2/)] 
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APPENDIX B. COMPONENTS OF MATRIX Q 


j la'p'p 


^[a>{f(a* + 2z>)-p>(3x>-2y>)) 

+R 2 [2x A + f{y 2 +z 2 ))^ 


2 + 2a‘] 

3 


~2^Y [“V +a 2 (2fz 2 - p 2 (2x 2 -ly 2 )) + 

R 2 (f(y 2 +z 2 )- 2 /)] 


6 







